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A  novel  stepped  flow  field  is  proposed  to  improve  the  performance  of  a  proton  exchange  membrane  fuel 
cell.  A  three-dimensional  model  that  accounts  for  activation  polarization,  ohmic  polarization  and  concen¬ 
tration  polarization  is  developed  to  achieve  this.  Numerical  results  demonstrate  that  overpotential  values 
decrease  in  the  order  t/act, c  >  0Po™m  >  >  ^ct,a  >  Oconc,c  >  ^conc,a  at  high  current  density.  Generally, 

the  concentration  polarization  can  be  ignored  because  of  its  low  value.  Using  a  stepped  flow  field  improves 
the  reactant  concentration  distribution,  local  current  density  distribution,  water  vapor  concentration  dis¬ 
tribution  and  cell  performance.  If  the  number  of  steps  approaches  infinity,  the  stepped  flow  field  extends 
to  a  tapered  field  that  has  the  lowest  cell  performance.  The  lower  the  number  of  steps  and  the  height  of 
the  outlet  channel,  the  greater  is  the  improvement  in  reactant  diffusion  to  the  porous  layers,  local  current 
density  distribution,  water  management  and  cell  performance. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  is  considered 
a  promising  alternative  power  source  to  the  internal  combustion 
engine  owing  to  its  high  efficiency,  competitive  power  density,  low- 
temperature  operation,  fast  start-up  and  zero  emission.  A  bipolar 
plate  is  a  vital  PEMFC  component  that  supplies  fuel  and  oxidant 
to  reactive  sites,  removes  reaction  products,  collects  the  current 
produced,  and  provides  mechanical  support  for  the  cell.  Bipolar 
plates  constitute  more  than  60%  of  the  weight  and  30%  of  the  total 
cost  of  a  fuel  cell  stack  [1,2]. 

To  obtain  high  PEMFC  performance,  various  flow  field  types  have 
been  proposed.  A  parallel  flow  field  is  a  conventional  type  of  flow 
field  for  which  the  channel  is  easily  fabricated,  but  the  cell  perfor¬ 
mance  is  relatively  low.  Hence,  some  new  flow  field  types  have  been 
investigated.  Hu  et  al.  [3]  investigated  the  effects  of  the  geometric 
parameters  of  a  parallel  flow  field  on  a  PEMFC  using  a  mathemati¬ 
cal  model.  They  found  that  fine  and  close  channels  were  beneficial 
for  cell  performance.  Yoon  et  al.  [4]  also  found  that  a  narrow  rib 
width  can  improve  cell  performance.  Relative  to  the  parallel  flow 
field,  a  serpentine  flow  field  yields  better  cell  performance  since 
the  gas  concentration  is  well  distributed  [5].  Weng  et  al.  [6]  com¬ 
pared  the  cell  performance  of  parallel  and  serpentine  flow  fields. 
Their  results  revealed  that  the  serpentine  flow  field  has  strong  con¬ 
vection  and  high  mass  transfer,  and  hence  better  cell  performance. 
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However,  serpentine  flow  fields  exhibit  large  pressure  drops  and 
water  flooding  can  occur  quite  easily  [1].  Wang  et  al.  [7]  analyzed 
the  effects  of  the  geometric  parameters  of  serpentine  flow  fields 
on  cell  performance.  Their  results  revealed  that  an  increase  in  the 
number  of  flow  channel  bends  can  improve  cell  performance,  that 
a  single  serpentine  flow  field  yields  better  performance  than  dou¬ 
ble  and  triple  serpentine  flow  fields,  and  that  cell  performance 
slowly  improves  with  increasing  flow  channel  width.  Compared 
to  parallel  and  serpentine  flow  fields,  an  interdigitated  flow  field 
yields  much  higher  cell  performance  since  reactant  transport  to 
the  diffusion  and  catalyst  layers  is  mainly  driven  by  forced  convec¬ 
tion  [6,8-10].  However,  interdigitated  flow  fields  suffer  from  very 
large  pressure  drops.  To  overcome  the  limitations  described,  some 
new  flow  field  types  have  been  investigated.  For  instance,  Kumar 
and  Reddy  [11]  and  Senn  and  Poulikakos  [12]  numerically  simu¬ 
lated  cell  performance  using  a  metal  foam  flow  field  and  observed 
a  more  uniform  local  current  density  distribution  and  higher  cell 
performance  compared  to  parallel  flow  fields.  Kuo  et  al.  [13]  numer¬ 
ically  investigated  PEMFC  performance  using  a  wave-like  flow  field 
and  observed  enhanced  transport  of  reactant  gases  through  the 
porous  layer,  increased  gas  flow  velocity,  and  improved  efficiency 
of  the  catalytic  reaction.  Consequently,  compared  to  a  parallel  flow 
field,  the  PEMFC  performance  for  a  wave-like  flow  field  was  notably 
higher.  Yan  and  co-workers  [14,15]  proposed  a  novel  baffle-blocked 
flow  field  design.  Baffles  are  arranged  in  the  cathode  flow  channel 
to  force  the  fuel  gas  to  flow  into  the  gas  diffusion  and  catalyst  lay¬ 
ers.  The  results  demonstrated  that  both  reactant  transport  and  cell 
performance  were  enhanced  by  the  presence  of  baffles  in  the  flow 
channel.  These  beneficial  effects  increased  with  increasing  width 
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Nomenclature 

A  area  (m)2 

As  specific  area  of  the  catalyst  layer  (m-1 ) 

C  molar  concentration  mol  (m-3 ) 

D  diffusion  coefficient  (m2  s-1 ) 

E  equilibrium  thermodynamic  potential  (V) 

F  Faraday  constant,  96487  C  (mol-1 ) 

H  height  (m) 

i  reaction  rate  (A  m-3 ) 

/  current  density  (A  m-2 ) 

K  permeability  (m2) 

L  length  (m) 

M  molar  mass  (kg  mol-1 ) 

n  electron  number  for  electrochemical  reactions 

N  number  of  steps 

p  pressure  (Pa) 

r  ratio  of  the  outlet  channel  height  to  inlet  channel 

height 

R  resistance  (£2  cm2);  gas  constant  (8.314 J mol-1  K-1) 

RH  relative  humidity 

S  source  term 

T  temperature  (K) 

u  velocity  vector  (ms-1) 

V  voltage  (V) 

w  velocity  (ms-1 ) 

W  width  (m) 

x,  y,  z  coordinate  direction  (m) 

Greek  symbols 

a  transfer  coefficient 

net  water  transfer  rate 
s  porosity 

g  stoichiometric  ratio 

r)  overpotential  (V) 

k  electrical  conductivity  (S  m-1 ) 

fi  fluid  viscosity  (kg  m-1  s-1 ) 

p  density  of  the  fluid  (kg  m-3 ) 

co  mass  fraction 

Subscripts  and  superscripts 
0  before  the  diffusion  layer 

a  anode 

act  activation 

av  average 

c  cathode 

cc  current  collector 

cone  concentration 

ct  catalyst  layer 

cell  fuel  cell 

ch  channel 

d  diffusion  layer 

eff  effective  value 

el  electron  transfer 

h  hydrogen 

in  inlet  channel 

pro  ion  transfer 

k  species 

L  limiting 

m  mass  and  membrane 

o  oxygen 

ohm  ohmic  polarization 

out  outlet  channel 


ref 

reference 

s 

specific 

u 

velocity 

w 

water 

and/or  number  of  baffles  in  the  tandem  array.  However,  there  is 
an  optimum  number  of  baffles  for  increasing  the  cell  performance. 
Based  on  a  fully  blocked  flow  field,  a  partially  blocked  flow  field 
that  can  obviously  improve  the  cell  performance  was  reported  by 
Soong  et  al.  [16]. 

The  flow  field  designs  described  above  can  significantly  improve 
the  cell  performance,  but  all  result  in  a  greater  pressure  drop  com¬ 
pared  to  a  parallel  flow  field.  Recently,  a  novel  tapered  flow  field 
was  proposed  by  Liu  et  al.  [17,18].  Their  numerical  results  indi¬ 
cate  that  this  tapered  flow  field  can  significantly  improve  the  cell 
performance  with  no  obvious  increase  in  pressure  drop. 

Here  we  report  on  a  novel  stepped  flow  field.  The  height  of  the 
flow  channel  decreases  in  a  stepwise  manner  along  the  main  flow 
direction.  If  the  number  of  steps  is  infinite,  the  stepped  flow  field 
extends  to  a  tapered  field.  The  objective  of  the  present  study  was  to 
numerically  investigate  PEMFC  performance  for  this  new  stepped 
flow  field  design.  To  carry  out  the  work,  a  three-dimensional,  single¬ 
phase  and  isothermal  numerical  model  accounting  for  activation 
polarization,  ohmic  polarization  and  concentration  polarization 
was  developed. 


2.  Model  description 


A  schematic  view  of  a  PEMFC  with  parallel  flow  fields  and  the 
computational  domain  is  shown  in  Fig.  1.  It  is  assumed  that  the 
PEMFC  structure  is  repeated  periodically  along  the  y-direction.  To 
save  computational  time,  a  typical  unit  shown  in  the  figure  is  taken 
as  the  computational  domain.  In  the  present  study,  stepped  flow 
fields  were  adopted  for  both  the  anode  and  cathode.  As  an  example, 
channels  with  three  steps  are  shown  in  Fig.  2.  In  fact,  the  number 
of  steps  is  a  factor  affecting  PEMFC  performance  that  is  analyzed  in 
the  following  sections. 

The  ratio  of  the  outlet  channel  height  to  the  inlet  channel  height, 
rch,  is  defined  as  follows: 


rch  = 


f^out 


(i) 


where  Hout  and  Hin  are  the  outlet  and  inlet  channel  height,  respec¬ 
tively. 


Computational  domain 


Current  collector 
Anode  diffusion  layer 
Anode  catalyst  layer 
Membrane 

Cathode  catalyst  layer 
Cathode  diffusion  layer 
Current  collector 


Gas  channel 


Fig.  1.  Schematic  view  of  a  PEMFC. 
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Fig.  2.  A  two-dimensional  PEMFC  in  the  x-z  plane  with  a  stepped  flow  field. 


Dry  air  is  fed  into  the  cathode  channel,  whereas  humidified 
hydrogen  is  supplied  to  the  anode  channel.  The  assumptions 
adopted  in  the  present  model  are  as  follows. 

(1)  The  fuel  cell  operates  under  steady-state  and  isothermal  con¬ 
ditions; 

(2)  The  gas  mixture  is  an  incompressible  ideal  fluid; 

(3)  The  gas  flow  in  the  channels  is  laminar; 

(4)  All  the  porous  zones  in  the  fuel  cell  domain  are  assumed  to  be 
isotropic  and  homogeneous,  and  the  membrane  is  considered 
impervious  to  reactant  gases;  and 

(5)  The  inlet  oxidant  is  not  humidified,  and  water  produced  on  the 
cathode  side  is  in  the  gaseous  state. 

2.1.  Governing  equations 

The  three-dimensional  model  consists  of  non-linear,  coupled 
partial  differential  equations  representing  the  conservation  of 
mass,  momentum  and  species.  The  equations  are  described  in  vec¬ 
tor  form  as  follows. 

The  mass  conservation  equation  can  be  expressed  as: 

V  •  (pu)  =  5m  (2) 

where  p  is  the  fluid  density,  u  is  the  velocity  vector  and  Sm  is 
a  source  term  for  electrochemical  reactions  in  the  catalyst  layers 
described  by: 

Anode  :  5m, a  —  5^  +  5w,a  (3) 

Cathode  :  5m, c  —  5o  +  5w,c  (4) 

The  source  terms  in  the  above  two  equations  will  be  given  in  the 
following  sections. 

The  momentum  conservation  equation  can  be  expressed  as: 

•  (puu)  =  -Vp+  •  (/xVu)  +  5u  (5) 

where  s  is  the  porosity  of  the  porous  medium,  p  is  the  pressure,  and 
5U  is  the  source  term,  which  is  as  follows  for  the  porous  medium: 


k,  and  Dk  eff  is  the  effective  diffusion  coefficient  of  species  /<,  which 
can  be  described  as  follows  in  the  porous  medium  [26]: 


(8) 

where  D/<  eff  is  the  reference  value  at  Tef f  and  Pref.  In  the  porous 
medium  region,  the  effective  diffusivity  of  the  species  can  be 
described  by  the  Bruggeman  model  [26]: 

Elk, eff  ~  5  '  E)k 

0) 

Sk  is  the  source  term  for  electrochemical  reactions  in  the  catalyst 
layers,  which  can  be  expressed  as: 

Sh  =  -jpMh 

(10) 

5w.a  = 

(ID 

S0  =  -^M0 

(12) 

S  ll+2ftM 

Jw,c  —  2  p  1V1W 

(13) 

The  cell  voltage  is  obtained  by  subtracting  all  overpotential  val¬ 
ues  from  the  equilibrium  thermodynamic  potential  as  follows  [19]: 


VCell  —  E  —  pact  —  Pohm  —  *1 cone  (14) 

where  E  is  the  equilibrium  thermodynamic  potential  that  can  be 
calculated  using  the  Nernst  equation  [20]: 

RT 

E  =  1.23  -  0.9  x  10 “3(T  -  298)  +  2.3  —  log(p2p0)  (15) 

Pact  is  the  activation  overpotential,  pohm  is  the  ohmic  overpoten¬ 
tial  and  pconc  is  the  concentration  overpotential.  These  parameters 
are  discussed  in  detail  as  follows. 

2.1.1.  Activation  overpotential 

The  activation  overpotential  comprises  anode  and  cathode  con¬ 
tributions.  The  relationship  between  the  overpotential  and  the 
reaction  rate  is  expressed  by  the  Butler-Volmer  equation: 


ic 


acncF 

Pact,c 


-  exp 


(1  —  ac)ncF 


RT 


Pact,c 


(17) 


2.1.2.  Ohmic  overpotential 

The  ohmic  overpotential  comprises  contributions  from  the 
resistance  to  electron  transfer  and  to  ion  transfer  [21].  Based  on 
Ohm’s  law,  it  can  be  expressed  as: 

^ohm  =  4m  +  Cm='(Rel+RPr°)  («) 


Su  =  -£eu  (6) 

The  species  conservation  equation  can  be  expressed  as: 

V  ■  (epuCfc)  =  V  ■  (pDfc  effVQ)  +  Sj<  (7) 

where  the  index  k  refers  to  different  species,  including  hydrogen, 
oxygen  and  water  vapor,  Ck  is  the  molar  concentration  of  species 


where  Rel  is  the  resistance  to  electron  transfer  and  Rpro  is  the  resis¬ 
tance  to  ion  transfer.  Rel  can  be  taken  as  a  constant,  but  is  generally 
difficult  to  predict  [21].  Based  on  the  computational  values  in  Ref. 
[22],  it  is  assumed  that  Rel  =  0.1  £2  cm2  in  the  present  model. 

The  resistance  to  ion  transfer  can  be  calculated  as: 

Rpr°  =^m  (19) 
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where  Hm  is  the  membrane  height  and  Km  is  the  membrane  phase 
conductivity. 


2A.3.  Concentration  overpotential 

The  concentration  overpotential  of  a  PEMFC  comprises  anode 
and  cathode  contributions,  which  can  be  calculated  as  [19]: 


where  Jl  is  the  limiting  current  density,  and  can  be  expressed  as 
[23]: 


_  nFDh^k,0 

Hd 


(21) 


Based  on  the  above  discussions,  the  cell  voltage  can  be  expressed 
as: 


VCell  —  E  ~  Pact, a  ~  Pact,c  ~  ^ohm  —  ^ohm  —  ^conc,a  —  Pconc,c  (22) 


2.2.  Boundary  conditions 


At  the  gas  channel  inlet,  the  species  concentrations  are  assumed 
to  be  constant.  The  inlet  velocities  are  given  by: 


^a,in  —  £a 


/ref  1  Am 
217  Qi,in  Ach 


^c,in  —  fc 


Ire  f  I  Am 
4F  Q),in  A ^ 


(23) 

(24) 


where  and  gc  are  the  stoichiometric  reactant  flow  ratios  for  the 
anode  and  cathode,  respectively,  and  are  defined  as  the  ratio  of  the 
reactant  quantity  supplied  to  the  reaction  quantity  generated  for  a 
specific  reference  current  density  Ire f.  Am  is  the  geometrical  area  of 
the  membrane  and  Ach  is  the  cross-sectional  area  of  the  gas  channel. 

A  local  one-way  assumption  is  adopted  to  provide  the  channel 
outlet  velocity  condition,  and  is  then  corrected  using  a  global  mass 
conservation  constraint  [24]. 

At  the  body  surface,  the  no-slip  condition  is  applied  for  the 
velocity  and  the  non-permeable  condition  is  applied  for  the  species 
mass  fraction. 


3.  Numerical  procedure 

The  governing  equations  and  the  boundary  conditions  are  dis¬ 
cretized  using  the  finite  volume  method.  The  SIMPLEC  algorithm 
[24]  is  used  for  coupling  the  velocity  and  pressure.  Since  all  gov¬ 
erning  equations  are  coupled  to  one  another,  they  should  be  solved 
simultaneously  using  an  iterative  method.  The  solution  is  consid¬ 
ered  to  be  convergent  when  the  relative  error  for  each  dependent 
variable  between  two  consecutive  iterations  is  less  than  1.0  x  10-5. 

The  grid  system  is  62  x  22  x  56.  To  simulate  local  transport  phe¬ 
nomena  in  the  fuel  cell,  the  grid  arrangement  in  the  z-direction 
is  non-uniform.  A  grid  independence  test  was  carried  out  for 
nine  grid  systems.  Cell  voltage  results  computed  using  the  model 
under  different  grid  systems  for  Jav  =  0.6  A  cm-2  are  summarized  in 
Table  1 .  Considering  both  accuracy  and  economics,  a  grid  system  of 
62  x  22  x  56  was  selected  for  the  present  study. 

The  parameters  used  in  the  present  model  are  listed  in  Table  2. 

4.  Results  and  discussion 

In  this  section,  p-Iav  curves  for  a  PEMFC  with  a  parallel  flow  field 
are  first  presented.  Then  the  performance  of  the  stepped  flow  field 
for  different  numbers  of  steps  is  compared.  Third,  the  effect  of  the 
outlet  channel  height  on  cell  performance  is  discussed.  Finally,  the 
pressure  drops  for  stepped  flow  fields  are  analyzed. 


Table  1 

Grid  independence  test  (Jav  =  0.6  A  cm-2). 


Stage 

Grid  size 

Vcell/V 

1 

52  x  22  x  56 

0.68360 

62  x  22  x  56 

0.68361 

72  x  22  x  56 

0.68362 

2 

62  x  12  x  56 

0.68347 

62  x  22  x  56 

0.68361 

62  x  34  x  56 

0.68370 

3 

62  x  22  x  40 

0.68355 

62  x  22  x  56 

0.68361 

62  x  22  x  64 

0.68362 

4 A.  Polarization  curve  and  overpotential  curves 

The  polarization  curve  and  overpotential  curves  for  a  PEMFC 
with  a  parallel  flow  field  are  shown  in  Figs.  3  and  4,  respectively. 
The  activation  polarization,  ohmic  polarization  and  concentration 
polarization  are  clearly  observed  in  Fig.  4  and  follow  the  general 
trends  reported  in  the  literature.  Overpotential  values  decrease  in 
the  order  pact,c  >  ^7^  >  ^hm  >  ?7act,a  >  Pconc,c  >  ^7conc,a  at  high 
current  density.  The  value  of  ^act,a  can  be  greater  than  that  of 
or  pfhm  at  low  current  density  owing  to  the  non-linear  relationship 
between  rjact}a  and  Jav.  In  comparison,  values  of  the  concentration 
overpotential  are  very  small.  For  instance,  when  /av  =  1.025  A  cm-2, 
Pconc,c  =  8.72  x  10-4  V  and  ^conc.a  =  2.63  x  10-3  V,  and  the  sum  of  the 


Table  2 

Parameters  used  in  the  model. 


Parameter 

Symbol 

Value 

Reference 

Gas  channel  length 

L 

0.06  m 

Gas  channel  width 

W 

1.0  x  10-3  m 

Gas  channel  height 

Hch 

2.0  x  10“3  m 

Diffusion  layer  height 

Hd 

2.54  x  10-4  m 

Catalyst  layer  height 

Hct 

2.87  x  10-5  m 

Membrane  height 

Hm 

2.3  x  10-4  m 

Land  area  width 

Wcc 

1.0  x  10-3  m 

Faraday’s  constant 

F 

96,487  C  mol-1 

Operation  temperature 

T 

353  K 

Anode/cathode  pressure 

Pa/Pc 

1/1  atm 

Electron  number  of  anode 

na 

4 

reaction 

Electron  number  of  cathode 

nc 

2 

reaction 

Fuel/air  stoichiometric  flow 

3/3 

[25] 

ratio 

Relative  humidity  of  inlet  fuel 

RHa 

100% 

[26] 

Relative  humidity  of  inlet  air 

RHc 

0 

[26] 

Oxygen  mass  fraction  of  inlet 

COo 

0.23 

air 

H2  diffusion  coefficient  at 

f^h.ref 

0.915  x  10“4  m2  s”1 

[25] 

reference  state 

02  diffusion  coefficient  at 

C^o.ref 

0.22  x  10“4  m2  s”1 

[25] 

reference  state 

Water  vapor  diffusion 

f^w,ref 

0.256  x  10-4  m2  s”1 

[25] 

coefficient  at  reference  state 

Anode  exchange  current 

^si'a.ref 

5.0  x  107  Am-3 

density  multiply  specific  area 

Cathode  exchange  current 

^s^c.ref 

120  Am-3 

density  multiply  specific  area 

Hydrogen  reference 

Qi.ref 

56.4  mol  m3 

[25] 

concentration 

Oxygen  reference 

Co, ref 

3.39  mol  m3 

[25] 

concentration 

Anode  transfer  coefficient 

aa 

0.5 

[26] 

Cathode  transfer  coefficient 

Ole 

0.5 

[26] 

Porosity  of  diffusion  layer 

Sd 

0.3 

[27] 

Porosity  of  catalyst  layer 

Set 

0.28 

[27] 

Absolute  permeability 

I< 

1.76  x  10-11  m2 

[27] 

Membrane  phase  conductivity 

/Cm 

17  Sm-1 

[25] 
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Fig.  3.  Polarization  curve  for  the  PEMFC.  FiS- 5-  Effect  of  N  on  the  cel1  Performance. 


two  values  is  approximately  4.4%  of  rjact<a.  Hence  the  concentration 
overpotentials  can  generally  be  ignored. 

In  particular,  the  anode  hydrogen  oxidation  reaction  (HOR)  is 
five  orders  of  magnitude  faster  than  the  cathode  oxygen  reduction 
reaction  (ORR),  and  therefore  one  may  consider  that  the  anode  acti¬ 
vation  overpotential  can  be  ignored.  However,  the  electron  number 
for  the  HOR  is  2,  which  is  half  of  the  number  for  the  ORR.  In  addi- 


Fig.  4.  Overpotential  vs.  average  current  density  for  (a)  higher  and  (b)  lower  over¬ 
potentials. 


tion,  the  model  equations  are  strongly  non-linear.  Hence,  the  anode 
activation  overpotential  is  not  negligible  in  terms  of  affecting  the 
cell  performance.  Similar  results  have  previously  been  reported 
[28]. 

4.2.  Effect  of  step  number  on  PEMFC  performance 

To  examine  the  effect  of  the  number  of  steps  on  PEMFC  perfor¬ 
mance,  I-V  curves  for  N=  1,  4,  and  12  at  rch  =  1/4  are  presented  in 
Fig.  5.  It  is  evident  that  the  effects  of  N  on  cell  performance  are  neg¬ 
ligible  at  low  current  density.  However,  variations  in  N  affect  cell 
performance  to  some  extent  at  high  current  density.  In  addition,  an 
increase  in  N  can  decrease  the  cell  performance.  Therefore,  we  can 
imagine  that  a  stepped  flow  field  will  extend  to  a  tapered  field,  as 
reported  in  Ref.  [17]  when  N  is  infinite,  which  exhibits  the  lowest 
cell  performance.  Hence,  a  lower  number  of  steps  is  beneficial  for 
cell  performance. 

To  explain  the  above  results,  distributions  of  the  average  oxygen 
concentration  in  the  cathode  catalyst  layer  and  the  cathode  current 
density  at  Jav  =  0.8  A  cm-2  are  shown  in  Figs.  6  and  7,  respectively. 
The  oxygen  concentration  and  cathode  current  density  gradually 
decrease  along  the  flow  direction  for  N=  12.  When  N=  1,  the  oxygen 


Fig.  6.  Distribution  of  oxygen  in  the  cathode  catalyst  layer  (Jav  =  0.8  A  cm-2 ):  (a)  N  =  1 
and  (b)  N=  12. 
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Fig.  7.  Distribution  of  the  cathode  current  density  (Jav  =  0.8  A  cm-2 ):  (a)  N  =  1  and  (b) 
N=  12. 


concentration  and  cathode  current  density  increase  significantly 
in  the  middle  of  the  channel  due  to  the  effect  of  the  step,  i.e.,  the 
reactant  gas  can  be  forced  more  effectively  into  the  porous  layers  for 
N=  1  compared  to  N=  12.  In  addition,  the  distributions  of  the  oxygen 
concentration  and  cathode  current  density  are  more  uniform  for 
N=1  than  for  N=  12. 

4.3.  Effect  of  outlet  channel  height  on  PEMFC  performance 

The  effects  of  rch  on  I-V  curves  are  shown  in  Fig.  8  for  N  =  3.  It 
is  evident  that  the  cell  performance  improved  with  decreasing  rch. 
As  in  the  discussion  above,  the  distribution  of  the  cathode  current 
density  is  used  to  illustrate  the  effects  of  rch  on  cell  performance. 
Fig.  9  shows  the  distribution  of  the  cathode  current  density  for 
rCh  =  4/16  and  13/16.  When  rch  is  smaller,  i.e.,  the  outlet  channel 
height  is  smaller,  steps  increase  the  local  current  density  and  lead 
to  an  improvement  in  cell  performance.  In  addition,  the  smaller  is 
rch,  the  more  uniform  is  the  current  density  distribution.  The  cell 
performance  therefore  improves. 

Water  management  is  very  important  for  PEMFCs.  To  illustrate 
the  high  performance  of  the  stepped  flow  field,  the  water  vapor 
distribution  at  the  cathode  is  discussed  in  detail.  Fig.  10  shows  the 


Fig.  8.  Effect  of  rch  on  the  cell  performance  (N  =  3). 


x  /m 


Fig.  9.  Distribution  of  the  cathode  current  density  (N= 3,  /av  =  0.8  A  cm-2):  (a) 
rch  =  4/16  and  (b)  rch  =  13/16. 


distribution  of  the  average  water  vapor  concentration  in  the  cath¬ 
ode  catalyst  layer  for  rch  =  4/16  and  13/16.  The  results  demonstrate 
that  when  rch  is  smaller,  i.e.,  the  outlet  channel  height  is  smaller, 
the  water  vapor  distribution  is  more  uniform  and  the  water  vapor 
concentration  is  lower.  These  phenomena  imply  that  a  stepped  flow 
field  with  small  outlet  channel  height  can  effectively  remove  water 
from  the  fuel  cell. 

4.4.  Pressure  drop  for  the  flow  field 

Finally,  we  consider  the  pressure  drop  for  flow  fields.  The  pres¬ 
sure  drop  between  the  cathode  inlet  channel  and  the  outlet  channel 
is  shown  in  Fig.  11.  As  illustrated  in  Fig.  11(a),  the  pressure  drop 
decreases  significantly  with  the  number  of  steps.  Fig.  11(b)  shows 
the  effect  of  the  channel  height  ratio,  rch,  on  the  pressure  drop.  A 
stepped  flow  field  with  a  small  height  increases  reactant  diffusion 
into  the  catalyst  layer  and  improves  water  removal,  but  increases 
the  pressure  drop.  Flowever,  the  maximum  pressure  drop  is  only 
21  Pa  in  the  typical  cases  considered  in  Fig.  11.  This  implies  a  small 
pressure  drop  for  use  of  a  stepped  flow  field. 


Fig.  10.  Distribution  of  the  average  water  vapor  concentration  in  the  cathode  cata¬ 
lyst  layer  (N= 3,  Jav  =  0.8  Acrrr2):  (a)  rch  =  4/16  and  (b)  rch  =  13/16. 
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Fig.  11.  Pressure  drop  for  the  stepped  flow  field:  effect  of  (a)  step  number  and  (b) 
outlet  channel  height. 


5.  Conclusions 

A  new  design  involving  a  stepped  flow  field  was  proposed.  The 
cell  performance  of  this  new  design  was  numerically  investigated. 
The  following  major  conclusions  can  be  drawn. 

(1)  The  new  model  developed  in  this  work  can  reasonably  sim¬ 
ulate  the  three  overpotentials.  Overpotential  decreases  in  the 
order  J)act,c  >  >  ^'hm  >  7?ac :t,a  >  Sconce  >  »?conc,a  at  high 

current  density.  Generally,  the  concentration  polarization  can 
be  ignored. 


(2)  The  stepped  flow  field  will  extend  to  a  tapered  field  if  the  num¬ 
ber  of  steps  approaches  infinity.  The  smaller  the  number  of 
steps,  the  more  favorable  is  the  improvement  in  reactant  dif¬ 
fusion  to  the  porous  layers,  local  current  density  distribution 
and  cell  performance. 

(3)  A  lower  outlet  channel  height  is  beneficial  in  improving  the  gas 
flow,  local  current  density  distribution,  water  management  and 
cell  performance. 

(4)  Use  of  a  stepped  flow  field  can  lead  to  a  significant  penalty  in 
terms  of  pressure  drop.  However,  the  pressure  drop  observed 
was  actually  low  and  the  cell  performance  was  high.  Hence,  this 
new  flow  field  can  be  widely  adopted  in  the  future. 

Finally,  it  should  be  noted  that  removal  of  liquid  water  from  the 

fuel  cell  is  very  important.  More  research  is  needed  to  analyze  the 

cell  performance  for  a  two-phase  and  non-isothermal  model. 
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